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The origin of ultra-high energy cosmic rays (UHECRs) 
and neutrinos is still a mystery. Hadronic acceleration the- 
ory suggests that they should originate in the same sources 
(astrophysical or cosmological), together with gamma- 
rays. While gamma-rays have been linked to astrophys- 
ical sources, no point source of UHECRs or neutrinos 
have been found so far. In this paper, the multimessen- 
ger combination of UHECRs and neutrinos as a new ap- 
proach to the high energy particle point source search is 
suggested. A statistical method for cross-correlation of 
UHECR and neutrino data sets is proposed. By obtain- 
ing the probability density function of number of neu- 
trino events within chosen angular distance from observed 
UHECRs, the number of neutrino events in the vicin- 
ity of observed UHECRs, necessary to claim a discov- 
ery with a chosen significance, can be calculated. Dif- 
ferent angular distances (bin sizes) are considered due to 
the unknown deflection of cosmic rays in galactic and in- 
tergalactic magnetic fields. Possible observed correlation 
of the arrival directions of UHECRs and neutrinos would 
provide a strong indication of hadronic acceleration the- 
ory. Correlation of both types of messengers with the 
location of certain sets of observed astrophysical objects 
would indicate sites of acceleration. Any systematic off- 
set in arrival directions between UHECRs and neutrinos 
may shed more light on magnetic field deflection of cos- 
mic rays. 



1 Introduction 

The origin of high energy particles coming from the Uni- 
verse is still not known. Currently, it is suggested that 
cosmic rays (protons and nuclei) up to about 10 15 eV 



gain their energy in acceleration by shocks from su- 
perno va explosions (see for example Na gano & Watson 
2000). For cosmic rays with energies around 10 16 17 eV, 
it has been suggested that acceleration can occur in the 
interaction of partic les with multiple supernova remnants 
dip & Axfordlll99lh . galacti c wind dVolk & Zirakashvili 
2004 . and microquasars dMirabelll2008l) . The highest 
energy charged cosmic rays are expected to be of extra- 
galacti c origin and ge t accelerated in jets of gamma-ray 
bursts (Waxman ll995blVietrj||l995h or active galactic nu- 
clei dBiermann & Strittmatterlll987t iRachen & Bierm ann 



19931) . Another option are so-called top-down scenar- 
ios: protons and neutrons, but also neutrinos and gamma- 
rays, are produced from quark and gluon fragmentati on of 
heavy exotic particles formed in the early Universe dHilll 
1983b ISchramm & Hilllll983l) . According to the theory 
of hadronic acceleration, the ultra-high energy cosmic ray 
(UHECR) flux is expected to be accompanied by associ- 
ated fluxes of gamma-rays and neutrinos from pion de- 
cays formed in the co llision of protons with photons (for 
example see IZasll2005l) . 

When a cosmic ray interacts with a particle in the 
Earth's atmosphere, a shower of particles is produced 
propagating towards the ground with almost the speed 
of light. These so-called "cosmic ray air showers" can 
be detected by arrays of Cherenkov detectors, fluorescent 
detectors or radio antennas. Low energy cosmic rays are 
constantly bombarding the Earth producing such show- 
ers, but cosmic rays of energies above 10 19 eV are very 
rare. They are expected to occur less than once a year on 
an area of one square kilometer and large shower arrays 
like the Pierre Auger Observatory (3000km 2 ) are needed. 
Moreover, cosmic rays are deflected by magnetic fields 
on their way from source to the Earth, thus their arrival 
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directions do not point back to their sources. The ex- 
pected deflection due to the Galactic magnetic field is 
in the order of few degrees for ultra-high energy protons 
dTinvakov & Tkachevl 120051) . Also, it is not yet known 
whether large deflections are caused by e xtra-galactic 



magn etic fields. According to some authors (Dol ag et al 



2005b they are in the order of one degree and bellow, 
but some other a nalysis are predicting far larger numbers 
( Sigl et af]|2004 ). The highest energy cosmic rays are ex- 
pected to be suppressed by interaction with the cosmic 
microwave background radiation, which i s known as the 
Greisen-Zatsepin-Kuzmin (GZK) cut-off dGreisen 1966; 
Zatsepin & Kuzmin 1966). Due to this effect, they can 
be observed only if they originate in sources closer than 
about 200Mpc. Note that the size of the GZK sphere 
is still under discussion, and significantly smaller GZK 
spheres are s uggested by some autho rs (see for example 
Beckerl2008l) . The HiRes experiment dAbbasi et alj2004 
measurements are in agreement with the predicted flux 
suppression of the highest energy co smic rays, however 
the results of the AGASA experiment dTakeda etal] 19981) 
show a n excess compared to the predicted GZK suppres- 
sion dBahcall & Waxmanl 120031) . The Pierre Auger col- 
laboration re ported recently observations consiste nt with 
GZK cut-off dThe Pierre Auger collaborafionll2008l) . 

Neutrinos are difficult to detect due to their small cross 
section for interaction with matter. They can interact 
within the Earth's atmosphere and trigger air showers. 
The probability for interaction increases with the path 
length, thus they are most likely to produce very inclined 
cosmic ray air showers. Those showers may be trig- 
gered anyw here along the travel ing path and also close to 



the ground (Ca pelle et al.lll998l) . However, a neutrino is 



much more likely to undergo charged current interaction 
with nuclei in the Earth's interior, producing a muon that 
emits Cherenkov radiation while moving through water 
or ice faster than the speed of light in the corresponding 
medium. This principle is used in current neutrino tele- 
scopes, such as IceCube and ANTARES. The IceCube 
neutrino telescope is currently being built at the South 
Pole. The final configuration should reach around 80 
strings with photomultiplier tubes (PMTs) in a 1km 3 ge- 
ometric volume. Currently, data are being taken with 40 
strings. ANTARES was recently finished in the Mediter- 
ranean Sea close to the French coast. Its final configura- 
tion consists of 12 strings with PMTs in 0.03km 3 of in- 




Figure 1: Aitoff skyplot of 27 UHECR events with 
E>57EeV (squares) observed by the Pierre Auger Obser- 
vatory and AGNs from the VCV catalog at D<75Mpc 
(diamonds) in equatorial coordinates. The color code 
represents the redshift of AGNs, with purple being 
the closest and red being the furthest objects. The 
dashed line represent the supergalacti c plane (plot after 
The Pierre Auger Collaboration! d2008l) ) 



strumented volume. The largest problem for observation 
of cosmic neutrinos is the large atmospheric flux of sec- 
ondary neutrinos coming from air showers. Comparing 
to UHECRs, the advantage of neutrinos as messengers is 
that they are not deflected in magnetic fields, so their ar- 
rival direction point back to their sources. The neutrino 
interaction cross section increases with increasing energy, 
so the highest energy neutrinos are more likely to interact 
with the detector, but also within the Earth's interior. This 
means that neutrinos above around IPeV can be observed 
only if coming from horizontal directions, so their path 
through the Earth interior is shorter. 



2 Cross-correlation of neutrinos 
and cosmic rays 

The Pierre Auger Observatory reported an anisotropy in 
the arrival directions of ultra high en ergy cosmic rays 
dThe Pierre Auger Collaborationll2008l) . Correlation with 
Active Galactic Nuclei (AG N) from the VCV catalog 
dVeron-Cettv & Veronl 120061) was the most significant 
for cosmic rays with energies higher than 57EeV and 
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Figure 2: Aitoff maps of combined sky visibility in equatorial coordinates. With the additional cut on zenith angles for 
UHECR observations, Z<60°, being taken into account, the field of view for the Pierre Auger Observatory is beneath 
declinations of 25° (on both panels: green). The left panel represents additionally the ANTARES neutrino telescope 
field of view (hatched). Visible declinations are beneath 47°. The right panel shows additionally the IceCube neutrino 
telescope field of view (hatched). Visible declinations (for up-going events) are above 0°. 



AGNs at distances less than 75Mpc. For this correla- 
tion, an angular separation between AGNs and cosmic 
rays larger than the point spread function of the detec- 
tor is considered, so magnetic deflection of cosmic rays 
is in this way partially taken into account. However, the 
mass composition of the observed events is not yet es- 
tablished, and the possible deflection in magnetic fields 
is not yet known. This means that in the case that the 
observed events are mostly nuclei, their arrival directions 
may significantly differ from the directions of their accel- 
eration sites. The anisotropy was estimated with a con- 
fidence level of 99%. The suggested correlation with 
local AGN sources mostly following the location of the 
supergalactic plane decreased in the following analysis 
( The Pierre Auger Collaborationl2009l) with 58 events ob- 
served at energies above 55EeV. The anisotropy of ob- 
served events remains with a confidence level of 99%. On 
the other side, no significant excess above the atmospheric 
neutrino flux was repor ted by the neutrino telescopes Ice- 
Cube and ANTARES (iceCube Collaboration: R. Abbasi 



20091: iToscano & the ANTARES collaborationll2009l) . 

Instead of searching for such a localized excess, neu- 
trino arrival direction can be correlated with the arrival di- 
rections of ultra high energy cosmic rays. In this paper, a 
method is suggested for such a correlation. The described 
analysis is based on the calculation of the probability den- 
sity function of the number of neutrino events within spe- 
cific angular distances from UHECRs. For a given set 



of observed cosmic ray arrival directions that are used as 
a catalog, Monte Carlo neutrino sets are generated, and 
angular distances between all cosmic rays and all MC 
neutrinos are calculated, (full sky correlation). In this 
way we obtain neutrino counts in chosen bins around the 
observed UHECRs, expected only from the atmospheric 
background neutrino flux. In other words, we calculate 
the probability that any neutrino count appears in the ob- 
served data. Usually, the bin size for this kind of analysis 
is the angular resolution of the detector, but due to the un- 
known deflection of cosmic rays in magnetic fields, we 
leave this parameter free. This allows our analysis to be 
sensitive for a potential systematic offset between arrival 
directions of ultra-high energy cosmic rays and neutrinos. 

An example for the calculation of the probability den- 
sity function of the number of neutrino events within spe- 
cific angular distances from UHECRs is presented. Fixed 
27 UHECRs events are used as catalog sources. Monte 
Carlo simulations with 1000 neutrinos on each generated 
sky map, randomly distributed over the full sky (decli- 
nations [90°,-90°] and right ascension [0°,360°]) are per- 
formed. Further, angular distances between all Monte 
Carlo neutrinos and 27 UHECRs are calculated, for each 
produced sky map. This provides the expected angu- 
lar distance distribution between neutrinos and cosmic 
rays in the absence of any correlation, over the whole 
sky. For a specific neutrino telescope, such as ANTARES 
or IceCube, instead of the whole sky, only declinations 
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Figure 3: Probability density function of the number of 
neutrino events on 1°,2°,5° and 10° angular distances 
from 27 fixed UHECR, for 1000 random background neu- 
trinos. 



visible by the chosen detector should be used, and in- 
stead of a completely random neutrino background dis- 
tribution, one that corresponds to the detector exposure 
should be applied. Also, the actually observed number 
of neutrino events should be used for Monte Carlo cal- 
culations. Correlation should be performed on a part of 
sky visible for both UHECR and neutrino detector. Fig. 
[2] shows the combined fields of view for the Auger ob- 
servatory and the ANTARES telescope (left panel) and 
the Auger observatory and the IceCube telescope (right 
panel). The Auger observatory is located at 35°S, and 
declinations up to 55° are observable. However, for the 
observed UHECR, an additional cut on zenith angle of 
60° is added, which means that declinations up to about 
25° are visible. The ANTARES telescope is located at 
43 °N, and for up-going events, declinations up to 47° are 
visible. The IceCube detector is at the South pole, so 
everything above 0° declination is in the field of view. 
IceCube is also investigating the possibility of extended 
sky search, by using also down-going events, and in that 
case all declinations above around -50° are observable 



(ILauer & the IceCube Collaborationll2009l) . 

Using this probability density function, the expected 
number of neutrino events in the vicinity of the observed 
ultra-high energy cosmic rays, necessary to claim a dis- 
covery with a chosen significance, can be calculated. 
Even though the expected angular resolution for IceCube 



and ANTARES is better than 1 degree, bin sizes of 1-10 
degrees are considered. This is because UHECRs may be 
deflected significantly in magnetic fields on their path to 
the Earth, thus the analysis is sensitive to a possible offset 
between arrival directions of cosmic rays and neutrinos. 

Probabilities for the expected number of background 
neutrinos, in bins around the 27 fixed UHECR events are 
shown in Fig. [3] For example, the probability to ob- 
serve 1 neutrino, within 27 bins of 1° is around 26%, and 
the probability to observe 5 neutrinos within those bins 
is around 3%. The probability to not have any neutrinos 
is around 14%, and to have more than 8 is close to zero. 
The same distributions are shown for bin sizes of 2,5 and 
10 degrees. This can be translated into discovery poten- 
tial and the number of neutrino events within a certain 
distance from observed cosmic rays, necessary to claim a 
discovery, can be obtained. 

Depending on Monte Carlo statistics, the number of 
neutrino events in the vicinity of cosmic rays, needed for 
a 3, 4 or 5<x discovery (correlation) claims can be estab- 
lished. This is shown in the left panel of Fig. |4]for 10 7 
Monte Carlo simulations, 27 fixed UHECRs, 1000 ran- 
dom neutrinos, and bin sizes of 1-10°. For example, cor- 
relation significance of 3cr is reached with about 1 1 neu- 
trino events within 27 bins of 1°, Act with 18 and 5<x with 
22. The right panel of Fig. |4]shows the percentage of neu- 
trino events necessary for 5<x discovery claim, depending 
on the bin size and the number of UHECRs. The bottom 
red line corresponds to magenta values presented in the 
left panel of Fig. |4] It is shown that in the case of 27 UHE- 
CRs, correlation significance of 5<x is reached with 2%- 
25% of neutrino events falling in 1-10° bins. However, 
with 100 UHECRs, it becomes more difficult to claim 5<x 
discovery for large bins, since 50-100% of neutrino events 
should be in 5-10° bins. This becomes even more extreme 
for 200 and 500 UHECRs, where more than 80% of neu- 
trino events should be already in 4° and 2° respectively, 
to claim a 5<x correlation. This stays the case even if a 
far larger amount of neutrinos is added (50000). How- 
ever, it does slightly depend on UHECRs sky distribution, 
since the exposure for neutrino detection depends on dec- 
lination. Anyhow, in the case of large number (>100) of 
observed UHECRs, it is suggested to first group those in 
clusters, and then apply the described method on centers 
of those UHECR clusters. 

Instead of the full (visible) sky, the UHECR-neutrino 
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Figure 4: Left panel: The number of neutrino events necessary to claim 3<x (red), 4<x (blue) and 5cr (magenta) 
correlation significance, as a function of the bin size, for the example of 27 fixed UHECRs and 1000 random neutrinos. 
Right panel: The percentage of neutrino events necessary to claim 5cr correlation significance, as a function of the bin 
size and number of UHECRs. 



cross correlation can also be applied only on certain re- 
gions. The most interesting regions, according to a cur- 
rent Auger results, are Centaurus A and the supergalactic 
plane due to the observed event excess. Also, the Virgo 
cluster from where no UHECR was yet observed, should 
be considered, since it might happen that due to a large 
magnetic field deflection, events are moved to other ar- 
rival directions. 

To test the method, fake neutrino sources that coincide 
with UHECR events are added to the Monte Carlo simu- 
lations. In this case, the expected number of neutrinos in- 
creases, proportional to the assumed number of neutrino 
sources coinciding with fixed UHECRs. This is shown in 
the left panel of Fig. [5] The luminosity of all fake sources 
is assumed to be 5 neutrinos (Gaussian-distributed around 
the given arrival direction of a cosmic ray), and the bin 
size is taken to be 1 degree. For example, the probability 
of having 5 neutrinos within 27 bins with one fake source 
is around 20%, compared to 5% with only background 
neutrinos. With 5 fake sources, the highest probability is 
to observe around 10 neutrinos in 27 bins. The probabil- 
ity to observe more than 20 neutrinos is beneath 1 %. With 
10 fake sources, around 20 neutrinos are most likely to be 
observed, and more than 30 are expected in less than 1 % 
of the cases. The right panel of Fig. [5]shows the expected 
distributions for 5 fake neutrino sources, with luminosi- 



ties 5, and exaggerated values of 10 and 20 neutrinos. If 
the added sources have higher luminosity, the expected 
number of neutrinos also increases. In other words, as 
expected, the method recognizes a stronger correlation, 
when the neutrino sources coincide with the observed cos- 
mic rays and when the neutrino sources are brighter. 

3 Discussion 

In the previous section, a statistical method for correla- 
tion of UHECR and neutrino data sets was described. 
Possible observed correlation of the arrival directions of 
those two types of messengers would provide a strong 
indication of hadronic acceleration theory. Correlation 
of both types of messengers with the location of cer- 
tain sets of observed astrophysical objects would indicate 
sites of acceleration. As we already mentioned, it is cur- 
rently believed that UHECR are originating in extragalac- 
tic sources, gamma-ray burst (GRBs) and active galactic 
nuclei (AGNs). While gamma-ray bursts (GRBs) are tran- 
sient extragalactic sources, and more suitable for correla- 
tion between gamma-rays and neutrinos (due to a signifi- 
cant delay in the arrival time between protons/nuclei and 
neutrinos), active galactic nuclei (AGNs), as extragalactic 
and steady sources, seem like a natural choice to explore 
the link between UHECRs and neutrinos. Assuming that 
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Figure 5: Left panel: Probability density function of the number of neutrino events within 1° angular distance from 
27 fixed UHECRs, for 1000 random neutrinos + fake neutrino sources (1,5 and 10) coinciding with UHECRs. Right 
panel: Probability density function of the number of neutrino events within 1° angular distance from 27 fixed UHE- 
CRs, for 1000 random neutrinos + 5 fake neutrino sources with luminosities 5, 10 and 20 neutrinos, coinciding with 
fixed UHECRs. 



there is a GZK cut-off, for an interpretation of correlation 
between arrival directions of UHECR and neutrinos, the 
limit should be set on AGNs within the GZK sphere, (left 
panel of Fig. [6}. However, we should keep on mind that 
AGNs with jets not pointing towards the Earth are pos- 
sible sources of UHECRs, since they can escape on the 
sides of jets, but neutrinos follow the direction of the jet 
and are possible observed only from blazar/BL Lac ob- 
jects, or ones so close by that they show blazar-like radi- 
ation, like it is the case with Centaurus A (right panel of 
Fig. 0. Due to UHECR excess observed in Cen A region 
by the Pierre Auger observatory, this region is interesting 
for UHECR/neutrino link investigatio n. Neutrino fluxes 
from C en A are recently calculated bv | Koers & Tinyakov 
d2008l) and iBecker & Biermannl d2009l) . but if neutrinos 
are indeed following the direction of jets, neutrino tele- 
scopes might not see the signal from this source. 

Any systematic offset in arrival directions between 
UHECRs and neutrinos may shed more light on mag- 
netic field deflection of cosmic rays. The Galactic mag- 
netic field has been shown to have both regular and ran- 
dom components. The regular component is believed to 
have a spiral s tructure and the corresponding global model 
was done by IStanevI ( 1997 ). According to this model, 
protons with an energy of about 4xl0 19 eV can be de- 
flected by the regular galactic magnetic field around 5°. 



Tinyakov & Tkachevl d2005l) calculated deflection in the 
random component to be about 0.2°- 1.5° d epending on 
the arr ival direction of a 4xl0 19 eV proton. iDolag et al 



(12005b used simulations of the large-scale structure to 
investigate magnetic fields in the intergalactic medium. 
They constructed full sky-maps of expected deflections 
for protons with energies of 10 20 eV and 4xl0 19 eV, and 
found that strong deflections are only present if UHE pro- 
tons cross a galaxy cluster, and even then stay beneath 
1°. On the other side. lSigl etail d2004 suggested that de- 
flection in extra-galactic magnetic fields can be of order 
20° up to 10 2() eV. They suggested that a region of strong 
magnetic field surrounding the observer would shield off 
UHECRs from sources outside of the magnetized region, 
so the observed flux would be dominated by a few closer 
sources and appear more anisotropic. 



4 Conclusions 

The hadronic acceleration theory predicts that ultra- 
high energy cosmic rays are accompanied by neu- 
trino and gamma-ray fluxes. While gamma-rays have 
been linked with astrophysical sources by observa- 
tions of telescopes like H.E.S.S. and MAGIC, and sev- 
eral attempts have been made to correlate them with 
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Figure 6: Left panel: Aitoff skymap in equatorial coordinates of 1205 AGNs from VCV catalog within 200Mpc 
(z<= 0.05) visible for the Pierre Auger Observatory. Colors are representing the redshift, with purple being the 
closest and red being the furthest objects. Right panel: Aitoff skymap in equatorial coordinates of BL Lac objects 
(confirmed, probable and questionable) + Cen A from VCV catalog within 200Mpc, visible for the Pierre Auger 
observatory. 



neutrinos (Hughev & the Ice Cube Collaboration 2007 
Ackermann etalJl2008t lBouwhuisll2008l) . the UHECR 



neutrino connection has not yet been examined. The rea- 
son is that both UHECRs and neutrinos have many down- 
sides as messengers. UHECRs are rare and they do not 
point back to their sources, since they are scrambled by 
galactic and intergalactic magnetic fields. Also, due to 
the interaction with the cosmic microwave background 
photons, they might be limited to the distance of about 
200Mpc or less. However, they are detectable with large 
shower arrays, like the Pierre Auger Observatory, which 
so far reported observation of few tens of events above 
55EeV. On the other side, cosmic neutrinos point back 
to their sources and their traveling distances are not lim- 
ited, but they interact rarely, which makes them difficult 
to detect. Also, they are for now lost in the huge atmo- 
spheric background of neutrinos coming from air show- 
ers, and so far none of the currently operating neutrino 
telescopes (IceCube, ANTARES) have observed neutrino 
excess above the atmospheric background flux. 

In this paper, the multimessenger combination of UHE- 
CRs and neutrinos as a new approach to the high en- 
ergy particle point source search is suggested. A statisti- 
cal method for cross-correlation of UHECR and neutrino 
data sets is proposed. By obtaining the probability density 
function of the number of neutrino events within specific 
angular distance from observed UHECRs, the number of 



neutrino events in the vicinity of observed ultra-high en- 
ergy cosmic rays, necessary to claim a discovery with a 
chosen significance, can be calculated. Different bin sizes 
are considered due to the unknown magnetic deflection 
of UHECRs. This multimessenger connection is crucial 
to indicate that hadronic acceleration actually happens in 
astrophysical sources, since only then neutrinos are to be 
observed coming from those sources. Through the cor- 
relation of arrival directions of those two very different 
messengers, we can gain additional insight into the non- 
thermal Universe. In the case of a correlation, we can 
learn more about the possible sources and physical pro- 
cesses that are responsible for the existence of ultra-high 
energy particles, such as hadronic acceleration of particles 
in sources. In the case of an offset between arrival direc- 
tions of cosmic rays and neutrinos, additional insight in 
UHECR deflection in magnetic fields can be gained. 
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